Clinical validation of genomic functional screen data: Analysis of observed BRCA1 variants in an unselected population cohort.

Functional assessment of genomic variants provides a promising approach to systematically examine the potential pathogenicity of variants independent of associated clinical data. However, making such conclusions requires validation with appropriate clinical findings. To this end, here, we use variant calls from exome data and BRCA1-related cancer diagnoses from electronic health records to demonstrate an association between published laboratory-based functional designations of BRCA1 variants and BRCA1-related cancer diagnoses in an unselected cohort of patient-participants. These findings validate and support further exploration of functional assay data to better understand the pathogenicity of rare variants. This information may be valuable in the context of healthy population genomic screening, where many rare, potentially pathogenic variants may not have sufficient associated clinical data to inform their interpretation directly.

Clinical outcomes of a genomic screening program for actionable genetic conditions.

PURPOSE: Three genetic conditions-hereditary breast and ovarian cancer syndrome, Lynch syndrome, and familial hypercholesterolemia-have tier 1 evidence for interventions that reduce morbidity and mortality, prompting proposals to screen unselected populations for these conditions. We examined the impact of genomic screening on risk management and early detection in an unselected population. METHODS: Observational study of electronic health records (EHR) among individuals in whom a pathogenic/likely pathogenic variant in a tier 1 gene was discovered through Geisinger's MyCode project. EHR of all eligible participants was evaluated for a prior genetic diagnosis and, among participants without such a diagnosis, relevant personal/family history, postdisclosure clinical diagnoses, and postdisclosure risk management. RESULTS: Eighty-seven percent of participants (305/351) did not have a prior genetic diagnosis of their tier 1 result. Of these, 65% had EHR evidence of relevant personal and/or family history of disease. Of 255 individuals eligible to have risk management, 70% (n = 179) had a recommended risk management procedure after results disclosure. Thirteen percent of participants (41/305) received a relevant clinical diagnosis after results disclosure. CONCLUSION: Genomic screening programs can identify previously unrecognized individuals at increased risk of cancer and heart disease and facilitate risk management and early cancer detection.

Patient-Centered Precision Health In A Learning Health Care System: Geisinger's Genomic Medicine Experience.

Health care delivery is increasingly influenced by the emerging concepts of precision health and the learning health care system. Although not synonymous with precision health, genomics is a key enabler of individualized care. Delivering patient-centered, genomics-informed care based on individual-level data in the current national landscape of health care delivery is a daunting challenge. Problems to overcome include data generation, analysis, storage, and transfer; knowledge management and representation for patients and providers at the point of care; process management; and outcomes definition, collection, and analysis. Development, testing, and implementation of a genomics-informed program requires multidisciplinary collaboration and building the concepts of precision health into a multilevel implementation framework. Using the principles of a learning health care system provides a promising solution. This article describes the implementation of population-based genomic medicine in an integrated learning health care system-a working example of a precision health program.

Linear assembly of a human centromere on the Y chromosome.

The human genome reference sequence remains incomplete owing to the challenge of assembling long tracts of near-identical tandem repeats in centromeres. We implemented a nanopore sequencing strategy to generate high-quality reads that span hundreds of kilobases of highly repetitive DNA in a human Y chromosome centromere. Combining these data with short-read variant validation, we assembled and characterized the centromeric region of a human Y chromosome.

Exome Sequencing-Based Screening for BRCA1/2 Expected Pathogenic Variants Among Adult Biobank Participants.

Importance: Detection of disease-associated variants in the BRCA1 and BRCA2 (BRCA1/2) genes allows for cancer prevention and early diagnosis in high-risk individuals. Objectives: To identify pathogenic and likely pathogenic (P/LP) BRCA1/2 variants in an unselected research cohort, and to characterize the features associated with P/LP variants. Design, Setting, and Participants: This is a cross-sectional study of adult volunteers (n = 50 726) who underwent exome sequencing at a single health care system (Geisinger Health System, Danville, Pennsylvania) from January 1, 2014, to March 1, 2016. Participants are part of the DiscovEHR cohort and were identified through the Geisinger MyCode Community Health Initiative. They consented to a research protocol that included sequencing and return of actionable test results. Clinical data from electronic health records and clinical visits were correlated with variants. Comparisons were made between those with (cases) and those without (controls) P/LP variants in BRCA1/2. Main Outcomes: Prevalence of P/LP BRCA1/2 variants in cohort, proportion of variant carriers not previously ascertained through clinical testing, and personal and family history of relevant cancers among BRCA1/2 variant carriers and noncarriers. Results: Of the 50 726 health system patients who underwent exome sequencing, 50 459 (99.5%) had no expected pathogenic BRCA1/2 variants and 267 (0.5%) were BRCA1/2 carriers. Of the 267 cases (148 [55.4%] were women and 119 [44.6%] were men with a mean [range] age of 58.9 [23-90] years), 183 (68.5%) received clinically confirmed results in their electronic health record. Among the 267 participants with P/LP BRCA1/2 variants, 219 (82.0%) had no prior clinical testing, 95 (35.6%) had BRCA1 variants, and 172 (64.4%) had BRCA2 variants. Syndromic cancer diagnoses were present in 11 (47.8%) of the 23 deceased BRCA1/2 carriers and in 56 (20.9%) of all 267 BRCA1/2 carriers. Among women, 31 (20.9%) of 148 variant carriers had a personal history of breast cancer, compared with 1554 (5.2%) of 29 880 noncarriers (odds ratio [OR], 5.95; 95% CI, 3.88-9.13; P < .001). Ovarian cancer history was present in 15 (10.1%) of 148 variant carriers and in 195 (0.6%) of 29 880 variant noncarriers (OR, 18.30; 95% CI, 10.48-31.4; P < .001). Among 89 BRCA1/2 carriers without prior testing but with comprehensive personal and family history data, 44 (49.4%) did not meet published guidelines for clinical testing. Conclusions and Relevance: This study found that compared with previous clinical care, exome sequencing-based screening identified 5 times as many individuals with P/LP BRCA1/2 variants. These findings suggest that genomic screening may identify BRCA1/2-associated cancer risk that might otherwise remain undetected within health care systems and may provide opportunities to reduce morbidity and mortality in patients.

Genomic characterization of large heterochromatic gaps in the human genome assembly.

The largest gaps in the human genome assembly correspond to multi-megabase heterochromatic regions composed primarily of two related families of tandem repeats, Human Satellites 2 and 3 (HSat2,3). The abundance of repetitive DNA in these regions challenges standard mapping and assembly algorithms, and as a result, the sequence composition and potential biological functions of these regions remain largely unexplored. Furthermore, existing genomic tools designed to predict consensus-based descriptions of repeat families cannot be readily applied to complex satellite repeats such as HSat2,3, which lack a consistent repeat unit reference sequence. Here we present an alignment-free method to characterize complex satellites using whole-genome shotgun read datasets. Utilizing this approach, we classify HSat2,3 sequences into fourteen subfamilies and predict their chromosomal distributions, resulting in a comprehensive satellite reference database to further enable genomic studies of heterochromatic regions. We also identify 1.3 Mb of non-repetitive sequence interspersed with HSat2,3 across 17 unmapped assembly scaffolds, including eight annotated gene predictions. Finally, we apply our satellite reference database to high-throughput sequence data from 396 males to estimate array size variation of the predominant HSat3 array on the Y chromosome, confirming that satellite array sizes can vary between individuals over an order of magnitude (7 to 98 Mb) and further demonstrating that array sizes are distributed differently within distinct Y haplogroups. In summary, we present a novel framework for generating initial reference databases for unassembled genomic regions enriched with complex satellite DNA, and we further demonstrate the utility of these reference databases for studying patterns of sequence variation within human populations.

An international effort towards developing standards for best practices in analysis, interpretation and reporting of clinical genome sequencing results in the CLARITY Challenge.

BACKGROUND: There is tremendous potential for genome sequencing to improve clinical diagnosis and care once it becomes routinely accessible, but this will require formalizing research methods into clinical best practices in the areas of sequence data generation, analysis, interpretation and reporting. The CLARITY Challenge was designed to spur convergence in methods for diagnosing genetic disease starting from clinical case history and genome sequencing data. DNA samples were obtained from three families with heritable genetic disorders and genomic sequence data were donated by sequencing platform vendors. The challenge was to analyze and interpret these data with the goals of identifying disease-causing variants and reporting the findings in a clinically useful format. Participating contestant groups were solicited broadly, and an independent panel of judges evaluated their performance. RESULTS: A total of 30 international groups were engaged. The entries reveal a general convergence of practices on most elements of the analysis and interpretation process. However, even given this commonality of approach, only two groups identified the consensus candidate variants in all disease cases, demonstrating a need for consistent fine-tuning of the generally accepted methods. There was greater diversity of the final clinical report content and in the patient consenting process, demonstrating that these areas require additional exploration and standardization. CONCLUSIONS: The CLARITY Challenge provides a comprehensive assessment of current practices for using genome sequencing to diagnose and report genetic diseases. There is remarkable convergence in bioinformatic techniques, but medical interpretation and reporting are areas that require further development by many groups.

Centromere reference models for human chromosomes X and Y satellite arrays.

The human genome sequence remains incomplete, with multimegabase-sized gaps representing the endogenous centromeres and other heterochromatic regions. Available sequence-based studies within these sites in the genome have demonstrated a role in centromere function and chromosome pairing, necessary to ensure proper chromosome segregation during cell division. A common genomic feature of these regions is the enrichment of long arrays of near-identical tandem repeats, known as satellite DNAs, which offer a limited number of variant sites to differentiate individual repeat copies across millions of bases. This substantial sequence homogeneity challenges available assembly strategies and, as a result, centromeric regions are omitted from ongoing genomic studies. To address this problem, we utilize monomer sequence and ordering information obtained from whole-genome shotgun reads to model two haploid human satellite arrays on chromosomes X and Y, resulting in an initial characterization of 3.83 Mb of centromeric DNA within an individual genome. To further expand the utility of each centromeric reference sequence model, we evaluate sites within the arrays for short-read mappability and chromosome specificity. Because satellite DNAs evolve in a concerted manner, we use these centromeric assemblies to assess the extent of sequence variation among 366 individuals from distinct human populations. We thus identify two satellite array variants in both X and Y centromeres, as determined by array length and sequence composition. This study provides an initial sequence characterization of a regional centromere and establishes a foundation to extend genomic characterization to these sites as well as to other repeat-rich regions within complex genomes.

Public knowledge of and attitudes toward genetics and genetic testing.

BACKGROUND: Variable health literacy and genetic knowledge may pose significant challenges to engaging the general public in personal genomics, specifically with respect to promoting risk comprehension and healthy behaviors. METHODS: We are conducting a multistage study of individual responses to genomic risk information for Type 2 diabetes mellitus. A total of 300 individuals were recruited from the general public in Durham, North Carolina: 60% self-identified as White; 70% female; and 65% have a college degree. As part of the baseline survey, we assessed genetic knowledge and attitudes toward genetic testing. RESULTS: Scores of factual knowledge of genetics ranged from 50% to 100% (average=84%), with significant differences in relation to racial groups, the education level, and age. Scores were significantly higher on questions pertaining to the inheritance and causes of disease (mean score 90%) compared to scientific questions (mean score 77.4%). Scores on the knowledge survey were significantly higher than scores from European populations. Participants' perceived knowledge of the social consequences of genetic testing was significantly lower than their perceived knowledge of the medical uses of testing. More than half agreed with the statement that testing may affect a person's ability to obtain health insurance (51.3%) and 16% were worried about the consequences of testing for chances of finding a job. CONCLUSIONS: Despite the relatively high educational status and genetic knowledge of the study population, we find an imbalance of knowledge between scientific and medical concepts related to genetics as well as between the medical applications and societal consequences of testing, suggesting that more effort is needed to present the benefits, risks, and limitations of genetic testing, particularly, at the social and personal levels, to ensure informed decision making.

Sequences associated with centromere competency in the human genome.

Centromeres, the sites of spindle attachment during mitosis and meiosis, are located in specific positions in the human genome, normally coincident with diverse subsets of alpha satellite DNA. While there is strong evidence supporting the association of some subfamilies of alpha satellite with centromere function, the basis for establishing whether a given alpha satellite sequence is or is not designated a functional centromere is unknown, and attempts to understand the role of particular sequence features in establishing centromere identity have been limited by the near identity and repetitive nature of satellite sequences. Utilizing a broadly applicable experimental approach to test sequence competency for centromere specification, we have carried out a genomic and epigenetic functional analysis of endogenous human centromere sequences available in the current human genome assembly. The data support a model in which functionally competent sequences confer an opportunity for centromere specification, integrating genomic and epigenetic signals and promoting the concept of context-dependent centromere inheritance.

Nonrandom X chromosome inactivation is influenced by multiple regions on the murine X chromosome.

During the development of female mammals, one of the two X chromosomes is inactivated, serving as a dosage-compensation mechanism to equalize the expression of X-linked genes in females and males. While the choice of which X chromosome to inactivate is normally random, X chromosome inactivation can be skewed in F1 hybrid mice, as determined by alleles at the X chromosome controlling element (Xce), a locus defined genetically by Cattanach over 40 years ago. Four Xce alleles have been defined in inbred mice in order of the tendency of the X chromosome to remain active: Xce(a) < Xce(b) < Xce(c) < Xce(d). While the identity of the Xce locus remains unknown, previous efforts to map sequences responsible for the Xce effect in hybrid mice have localized the Xce to candidate regions that overlap the X chromosome inactivation center (Xic), which includes the Xist and Tsix genes. Here, we have intercrossed 129S1/SvImJ, which carries the Xce(a) allele, and Mus musculus castaneus EiJ, which carries the Xce(c) allele, to generate recombinant lines with single or double recombinant breakpoints near or within the Xce candidate region. In female progeny of 129S1/SvImJ females mated to recombinant males, we have measured the X chromosome inactivation ratio using allele-specific expression assays of genes on the X chromosome. We have identified regions, both proximal and distal to Xist/Tsix, that contribute to the choice of which X chromosome to inactivate, indicating that multiple elements on the X chromosome contribute to the Xce.

Evidence for sequence biases associated with patterns of histone methylation.

BACKGROUND: Combinations of histone variants and modifications, conceptually representing a histone code, have been proposed to play a significant role in gene regulation and developmental processes in complex organisms. While various mechanisms have been implicated in establishing and maintaining epigenetic patterns at specific locations in the genome, they are generally believed to be independent of primary DNA sequence on a more global scale. RESULTS: To address this systematically in the case of the human genome, we have analyzed primary DNA sequences underlying patterns of 19 different methylated histones in human primary T-cells and patterns of three methylated histones across additional human cell lines. We report strong sequence biases associated with most of these histone marks genome-wide in each cell type. Furthermore, the sequence characteristics for such association are distinct for different groups of histone marks. CONCLUSIONS: These findings provide evidence of an influence of genomic sequence on patterns of histone modification associated with gene expression and chromatin programming, and they suggest that the mechanisms responsible for global histone modifications may interpret genomic sequence in various ways.

Composition and organization of active centromere sequences in complex genomes.

BACKGROUND: Centromeres are sites of chromosomal spindle attachment during mitosis and meiosis. While the sequence basis for centromere identity remains a subject of considerable debate, one approach is to examine the genomic organization at these active sites that are correlated with epigenetic marks of centromere function. RESULTS: We have developed an approach to characterize both satellite and non-satellite centromeric sequences that are missing from current assemblies in complex genomes, using the dog genome as an example. Combining this genomic reference with an epigenetic dataset corresponding to sequences associated with the histone H3 variant centromere protein A (CENP-A), we identify active satellite sequence domains that appear to be both functionally and spatially distinct within the overall definition of satellite families. CONCLUSIONS: These findings establish a genomic and epigenetic foundation for exploring the functional role of centromeric sequences in the previously sequenced dog genome and provide a model for similar studies within the context of less-characterized genomes.

Effects of sequence variation on differential allelic transcription factor occupancy and gene expression.

A complex interplay between transcription factors (TFs) and the genome regulates transcription. However, connecting variation in genome sequence with variation in TF binding and gene expression is challenging due to environmental differences between individuals and cell types. To address this problem, we measured genome-wide differential allelic occupancy of 24 TFs and EP300 in a human lymphoblastoid cell line GM12878. Overall, 5% of human TF binding sites have an allelic imbalance in occupancy. At many sites, TFs clustered in TF-binding hubs on the same homolog in especially open chromatin. While genetic variation in core TF binding motifs generally resulted in large allelic differences in TF occupancy, most allelic differences in occupancy were subtle and associated with disruption of weak or noncanonical motifs. We also measured genome-wide differential allelic expression of genes with and without heterozygous exonic variants in the same cells. We found that genes with differential allelic expression were overall less expressed both in GM12878 cells and in unrelated human cell lines. Comparing TF occupancy with expression, we found strong association between allelic occupancy and expression within 100 bp of transcription start sites (TSSs), and weak association up to 100 kb from TSSs. Sites of differential allelic occupancy were significantly enriched for variants associated with disease, particularly autoimmune disease, suggesting that allelic differences in TF occupancy give functional insights into intergenic variants associated with disease. Our results have the potential to increase the power and interpretability of association studies by targeting functional intergenic variants in addition to protein coding sequences.

Uncoupling of genomic and epigenetic signals in the maintenance and inheritance of heterochromatin domains in fission yeast.

Many essential aspects of genome function, including gene expression and chromosome segregation, are mediated throughout development and differentiation by changes in the chromatin state. Along with genomic signals encoded in the DNA, epigenetic processes regulate heritable gene expression patterns. Genomic signals such as enhancers, silencers, and repetitive DNA, while required for the establishment of alternative chromatin states, have an unclear role in epigenetic processes that underlie the persistence of chromatin states throughout development. Here, we demonstrate in fission yeast that the maintenance and inheritance of ectopic heterochromatin domains are independent of the genomic sequences necessary for their de novo establishment. We find that both structural heterochromatin and gene silencing can be stably maintained over an ~10-kb domain for up to hundreds of cell divisions in the absence of genomic sequences required for heterochromatin establishment, demonstrating the long-term persistence and stability of this chromatin state. The de novo heterochromatin, despite the absence of nucleation sequences, is also stably inherited through meiosis. Together, these studies provide evidence for chromatin-dependent, epigenetic control of gene silencing that is heritable, stable, and self-sustaining, even in the absence of the originating genomic signals.

Analysis of DNA methylation in a three-generation family reveals widespread genetic influence on epigenetic regulation.

The methylation of cytosines in CpG dinucleotides is essential for cellular differentiation and the progression of many cancers, and it plays an important role in gametic imprinting. To assess variation and inheritance of genome-wide patterns of DNA methylation simultaneously in humans, we applied reduced representation bisulfite sequencing (RRBS) to somatic DNA from six members of a three-generation family. We observed that 8.1% of heterozygous SNPs are associated with differential methylation in cis, which provides a robust signature for Mendelian transmission and relatedness. The vast majority of differential methylation between homologous chromosomes (>92%) occurs on a particular haplotype as opposed to being associated with the gender of the parent of origin, indicating that genotype affects DNA methylation of far more loci than does gametic imprinting. We found that 75% of genotype-dependent differential methylation events in the family are also seen in unrelated individuals and that overall genotype can explain 80% of the variation in DNA methylation. These events are under-represented in CpG islands, enriched in intergenic regions, and located in regions of low evolutionary conservation. Even though they are generally not in functionally constrained regions, 22% (twice as many as expected by chance) of genes harboring genotype-dependent DNA methylation exhibited allele-specific gene expression as measured by RNA-seq of a lymphoblastoid cell line, indicating that some of these events are associated with gene expression differences. Overall, our results demonstrate that the influence of genotype on patterns of DNA methylation is widespread in the genome and greatly exceeds the influence of imprinting on genome-wide methylation patterns.

Organization and molecular evolution of CENP-A--associated satellite DNA families in a basal primate genome.

Centromeric regions in many complex eukaryotic species contain highly repetitive satellite DNAs. Despite the diversity of centromeric DNA sequences among species, the functional centromeres in all species studied to date are marked by CENP-A, a centromere-specific histone H3 variant. Although it is well established that families of multimeric higher-order alpha satellite are conserved at the centromeres of human and great ape chromosomes and that diverged monomeric alpha satellite is found in old and new world monkey genomes, little is known about the organization, function, and evolution of centromeric sequences in more distant primates, including lemurs. Aye-Aye (Daubentonia madagascariensis) is a basal primate and is located at a key position in the evolutionary tree to study centromeric satellite transitions in primate genomes. Using the approach of chromatin immunoprecipitation with antibodies directed to CENP-A, we have identified two satellite families, Daubentonia madagascariensis Aye-Aye 1 (DMA1) and Daubentonia madagascariensis Aye-Aye 2 (DMA2), related to each other but unrelated in sequence to alpha satellite or any other previously described primate or mammalian satellite DNA families. Here, we describe the initial genomic and phylogenetic organization of DMA1 and DMA2 and present evidence of higher-order repeats in Aye-Aye centromeric domains, providing an opportunity to study the emergence of chromosome-specific modes of satellite DNA evolution in primate genomes.

Allele-specific distribution of RNA polymerase II on female X chromosomes.

While the distribution of RNA polymerase II (PolII) in a variety of complex genomes is correlated with gene expression, the presence of PolII at a gene does not necessarily indicate active expression. Various patterns of PolII binding have been described genome wide; however, whether or not PolII binds at transcriptionally inactive sites remains uncertain. The two X chromosomes in female cells in mammals present an opportunity to examine each of the two alleles of a given locus in both active and inactive states, depending on which X chromosome is silenced by X chromosome inactivation. Here, we investigated PolII occupancy and expression of the associated genes across the active (Xa) and inactive (Xi) X chromosomes in human female cells to elucidate the relationship of gene expression and PolII binding. We find that, while PolII in the pseudoautosomal region occupies both chromosomes at similar levels, it is significantly biased toward the Xa throughout the rest of the chromosome. The general paucity of PolII on the Xi notwithstanding, detectable (albeit significantly reduced) binding can be observed, especially on the evolutionarily younger short arm of the X. PolII levels at genes that escape inactivation correlate with the levels of their expression; however, additional PolII sites can be found at apparently silenced regions, suggesting the possibility of a subset of genes on the Xi that are poised for expression. Consistent with this hypothesis, we show that a high proportion of genes associated with PolII-accessible sites, while silenced in GM12878, are expressed in other female cell lines.

Comparative analysis of the primate X-inactivation center region and reconstruction of the ancestral primate XIST locus.

Here we provide a detailed comparative analysis across the candidate X-Inactivation Center (XIC) region and the XIST locus in the genomes of six primates and three mammalian outgroup species. Since lemurs and other strepsirrhine primates represent the sister lineage to all other primates, this analysis focuses on lemurs to reconstruct the ancestral primate sequences and to gain insight into the evolution of this region and the genes within it. This comparative evolutionary genomics approach reveals significant expansion in genomic size across the XIC region in higher primates, with minimal size alterations across the XIST locus itself. Reconstructed primate ancestral XIC sequences show that the most dramatic changes during the past 80 million years occurred between the ancestral primate and the lineage leading to Old World monkeys. In contrast, the XIST locus compared between human and the primate ancestor does not indicate any dramatic changes to exons or XIST-specific repeats; rather, evolution of this locus reflects small incremental changes in overall sequence identity and short repeat insertions. While this comparative analysis reinforces that the region around XIST has been subject to significant genomic change, even among primates, our data suggest that evolution of the XIST sequences themselves represents only small lineage-specific changes across the past 80 million years.